Pressure and intermittency in passive vector turbulence 
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We investigate the scaling properties a model of passive vector turbulence with pressure and in 
the presence of a large-scale anisotropy. The leading scaling exponents of the structure functions are 
proven to be anomalous. The anisotropic exponents are organized in hierarchical families growing 
without bound with the degree of anisotropy. Nonlocality produces poles in the inertial-range 
dynamics corresponding to the dimensional scaling solution. The increase with the Peclet number 
of hyperskewness and higher odd-dimensional ratios signals the persistence of anisotropy effects also 
in the inertial range. 
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In the last five years, much progress has been achieved in understanding of the physical origin of intermittency 
and anomalous scaling in fluid turbulence. The study of the passive advection of a scalar quantity (temperature, 
concentration of an impurity etc) by a random Gaussian field, white in time and self-similar in space, the rapid-change 
model jjj, has played a crucial role in these developments. There, it was for the first time possible to demonstrate 
the emergence of anomalous scaling and to relate it to the existence of statistical conservation laws of the dynamics 
v./ [ |^|,D . Such conservation laws appear as scaling zero modes of the hierarchy of linear operators governing the inertial- 
range dynamics of the correlation functions. The dominance of such zero modes has been confirmed in numerical 
experiments and demonstrated analytically by means of the renormalization group and operator product expansion 
||[3) . The importance of these results is seen from the fact that statistical invariants of dynamics are features of the 
passive advection by more general classes of turbulent velocity fields |jj . A review and more references can be found 
t-h ; in Ref. @. 

^ , The first-principle model of fluid turbulence is the Navier-Stokes equation, a nonlinear integro-differential equation 
for the solenoidal vector velocity field. Whether the relation between statistical invariants and anomalous scaling can 
be generalized to the nonlinear dynamics of the advecting velocity field remains an open and challenging issue. Setting 
^v^j ' aside the problem of the nonlinearity of interactions, the scope of the present Letter is to investigate how nonlocality 
\ affects scaling. 

t-H ■ In the Navier-Stokes equation nonlocality arises through the pressure term fixed by the incompressibility condition 

for the velocity field. Here, we consider the passive advection of an incompressible vector field u according to the 

£^ ■ most general dynamics consistent with the Galilean invariance: 
• 1— ( ' 

^ ! 9 t u + (v9)u-i(u^)v + 9P = /tAu + f, (1) 

. I-h I where V is the pressure and k is the diffusivity coefficient. The advection field v is specified by the Kraichnan ensemble. 

It models a homogeneous and isotropic nonintermittent velocity field and, due to Gaussianity, is determined completely 
by the second-order structure function 

S a! 3(v) = («„(x, t) Mx,i') - v/siy,?))) =D5(t- f)M £ [(d - 1 + - £f a f p ] (2) 

with r = x — y and f = r/\r\. Here < £ < 2 is a kind of Holder exponent which measures the "roughness" of the 
velocity field. In the renormalization group approach, it plays the same role as the parameter e = 4 — d in the theory 
of critical phenomena. 

The vector field u is stirred by an incompressible forcing f varying only over large spatial scales. We assume the 
forcing f to be delta-correlated in time and Gaussian in space with the covariance 

(f a ( X )f p (y))=5(t-t')F aP (r/L). (3) 

The parameter L is the integral scale of the stirring, and F a p is a function which rapidly decays as r/L — > 00 and 
tends to a constant as r/L — > 0. 

As in the Navier-Stokes equation, the pressure V is fixed by requiring Eq. ([!]) to be consistent with the incompress- 
ibility condition: 
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AV(x,t) = (A-l)(d a Vf,)(df,U a ). (4) 

The parameter A modules the degree of nonlocality in the problem. The case A = 1, known as the Kazantsev- 
Kraichnan kinematic dynamo model, represents the limit of local interactions. It describes the early stage of the 
growth of a magnetic seed in a turbulent fluid. However, the exponential growth of the magnetic field occurs only 
if £ is large enough, in particular for £ > I in three dimensions |J. For smaller 4, a steady state sets in exhibiting 
anomalous scaling behavior already for the pair correlation function. The corresponding exponents have been found 
exactly in the isotropic H and anisotropic [p~0| p~2|| sectors. 

As A deviates from unity, the solution of the Poisson equation (Q) introduces in the dynamics a nonlocal integral 
term. Cases of special interest are A — and A = — 1. For A — the stretching term (u • 9)v vanishes and 
the dynamics becomes invariant under translations u — > u + const. The kinetic energy (u 2 (x, t)) is conserved, as 
in the passive scalar and stochastic Navier-Stokes equation; see Ref. The case A = — 1 corresponds to the 

linearization of the Navier-Stokes equation around the rapid-change velocity field; in two dimensions it was studied 
in Ref. @. 

We are interested in the scaling properties of the correlation functions of the field u in the inertial range, L >> 
|r| >> t = (k/D) 1 '^, where i is the dissipative scale. The equal-time pair correlation function is governed by a 
linear mtegro-differential equation, 

d t C a p = T>1 S pCys, C a p(r,t) = {u a (x,t)up(y,t)) , (5) 

where the summation over repeated indices is implied. The explicit form of the inertial operator 2^a, of scaling 
dimension — 2 + 4, can be derived using functional integrations by parts JjiJ or the Dyson- Wyld equation pju]]. It is 
not given here for the sake of brevity. Under the assumption that a steady state sets in, the inertial-range isotropic 
solution is sought in the time-independent scaling form: 

C a0 (r) =C\t\C* [{d-l + ^)5 aP -C, 2 f a rp\. (6) 

After the adoption of the Anstaz (^), we define the integrals in Eq. (|^) in the spirit dimensional regularization. Taking 
the trace then gives the desired equation for £ 2 : 



VHC 1S = DC |r|- 2 +^ {d ^Z+f 2) C2,A)= 0, 

1(2 + 4 - l ) 



(7) 



where 
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T(C 2 , £,A) = -(A- l) 2 (d + 1)4 C2 + (C 2 + £ - 2) [d 2 (C 2 + A 2 4) + d (-C 2 + Cf + A 2 i{t- 1)) - (C 2 -AO 

iU 1)(Af a d n r(i + ^/2)r(i + d/2 + C2/2)r(2-e/2-C2/2)r(i + d/2) 
[A )(Ak A ' r(d/2 + 4/2 + C2/2)r(-C2/2)r(2 + d/2-4/2) ' w 

Equation (^) has always a simple root for £2 = 2 — d — 4 and for all A ^ 0,1a simple pole for (2 — 2 — 4- It is worth 
stressing that the pole arises because of the nonlocal terms in the inertial operator. The scaling | r | 2 ^ is marginal and 
produces logarithmic divergences at the edges of the inertial range which are reflected in the presence of the pole. In 
order to understand its physical meaning it is instructive to compare Eq. (Q) with the passive scalar case [^|J|]. The 
inertial operator governing the pair correlation function again has dimension —2 + 4; its action onto a power function 
yields 

p| r |-2+ e+C2 = D C | r |-2+ 5+C 2 (rf _ 1) C2 ( C2 + ^ + d _ 2 ). (9) 

The root (2 — 2 — d — 4 is again present. It corresponds to the dimensional solution of the equation with the right-hand 
side is proportional to 8(r) (or a constant in momentum space). The corresponding scaling behavior is admissible 
at scales much larger than the integral one The second root C2 = 0, a constant zero mode, is a consequence of 
energy conservation and as such it is also a root of Eq. (0) for ,4 = 0. Finally, setting £ 2 = 2 — 4 in Eq. (||) solves 
the equation with a constant right-hand side (or £(k) in the momentum space). The pole for £ 2 — 2 — 4 m Eq. (|^) 
thus corresponds to the scaling solution obtained by matching the forcing at large scales. Its cancellation for A = 
indicates the existence of a constant flux solution in the presence of kinetic energy conservation. 

The transcendental equation T(£ 2 , 4, A) — has infinitely many solutions. The leading admissible exponent behaves 
as C2 — O(0 for small 4- It can be obtained from Eq. (||) within perturbation theories in 1/d or 4: 
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>+_ A 2c,A 2 Z \ ( A - !) r (± + ^/ 2 ) r I 1 + C_? - €M + l)M 3 g - .4 2 g + -4g + -A - g + 1) 1 , nn , , 2 , 

C 2 - ^ + — j _____ _____ |+^i/dj 

= V+VhK - M^+^-S)* ^ + I) 2 + (* - « + ^ + 2-4 + 3)} + 0(e). (10) 

The nonlocality parameter A enters the perturbative expansions from the first order. The fact reflects the balance 
between local and nonlocal contributions in the zero-mode equation. It is worth noting that up to first order in 1/d 
the exponent of the linearized Navier-Stokes and of the magnetic, local, model coincide. 

The next-to-leading correction exponent has the form = 2 - £ + A(A - l)(d + 2)£/(d 2 + A 2 + Ad - 3) + 0(£ 2 ). 

Nonperturbative solutions of Eq. ([?]) can be obtained analytically only in the case A = 0, where = and 
(2 = 2 — £ exactly, and in the local case .4 = 1 when it reduces to a third-order algebraic equation. Besides the 
subleading solution = 2 — £, the result of Ref. ||] is recovered. In general, the behaviour of the exponents 
depends qualitatively on the value of A. Below we focus our attention only on the behaviors of the leading exponent 
Q2 vs £ for d — 3 and in the physically interesting case A E [—1,1]. More detailed analysis will be reported elsewhere. 

When —0.581 < A < 0.613, the solution £2" exists and is strictly negative within the entire interval < £ < 2, 
except for the case A — when Q' vanishes identically. For the other values of A, the real solution exists only in the 
subinterval < £ < £ c , (£ c (-4, d) < 2 is some critical value) where it is again strictly negative. For £ > £ c , the solution 
(2 coalesces with the closest unphysical branch, (2 = —d + 0(£), an d they both become complex: the effect known 
for the magnetic model where £ c (l, 3) = 1. It was argued in Refs. J|,[_l that this critical value of £ coincides with 
the threshold for the dynamo effect (exponential growth of the pair correlation function). 

It is worth noticing that in the small subinterval, —0.613 < A < —0.581, both solutions again become real near to 
£ = 2 (see Fig. 1). 

Coefficients of the expansions in £ for diverge at (d 2 + A 2 + Ad — 3) = (this can happen only for d < 2, 
and for d = 2 only for A = —1). Nonperturbative analysis shows that in the region [d 2 + A 2 + Ad — 3) < 0, the 
effective diffusivity coefficient is negative at large scales, which makes the system unstable with respect to any small 
perturbation. In the limit (d 2 + A 2 + Ad — 3) — > 0, the critical value £ c vanishes, and the solution for the leading 
inertial-range exponent disappears. We thus conclude that the solution (2 = — £, reported in Ref. |lj| for d = 2 and 
A = —1, in fact corresponds to the large-scale exponent C2 = 2 — d — £ mentioned above and not to the inertial range. 

Nonperturbative analysis can be repeated for the anisotropic sectors of the pair correlation function; see Refs. 
]To| |l2"t for A = 1 and [l||,|l4| for .4 = 0. Anisotropic contributions can be excited, e.g., in the presence of a constant 
background field jl0],[ll | . The correlation can be expanded in components of definite total angular momentum j with 
radial dependence of the powerlike form (parity considerations allow only even j in this expansion) : 

00 

C a0 (r)^J2C a0 (j) \r\^ Pjiv). (11) 

However, up to first order in £ these exponents can be most conveniently calculated using the renormalization group 
and operator product expansion. We refer the reader to the literature ]|,[|[ll],[l7| for details and present only the 
result: 
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(A 2 + d 2 + Ad- 3)(2j + d - 6)(2j + d- 4)(2j + d - 2) 
1 Afd + 2j 2 d 2 - 10j 3 - 22j 2 rf - 6jd 2 + 35j 2 + 38jd + 6d 2 - 50j - 24d + 24).4 2 

+(j - 2)(j 2 + 2jd + d 2 - 5j - 4d+4)[2A(j +d- 3) + (j - 3)]}, (12) 



where j > 2. The dependence on the angular momentum indicates the existence of an hierarchy of exponents growing 
with j without an upper bound. The result is in agreement with previous findings for the passive advection of the 



scalar |17|] and magnetic ]10|-|l2| fields and with the nonlocal case .4 = J13 14|. 

The renormalization group and operator product analysis is also extended to correlation functions of arbitrary 
order. The key role is played by the dimensions A[n,j] associated with the jth rank tensor composite operators 
u ai . . . u a . (u 2 )', where n = 21 + j. Let T> np (x, y) be some equal-time scalar correlation function of the n th and pth 
powers of the field u. For the leading term in the j th anisotropic sector we obtain: 

2W*,y) « (\r\/e)- A ^-^, jp (| r | /jL )A„ 3 p. ( f )) (13) 
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where jk — or 1 for even or odd k, respectively; cf. [p|Jl 7|] for the scalar and for the magnetic cases. Up to order 
0(£), we have obtained 

A[n ' J] = T 2(^+^ + ^-3) + ° K } - (M) 

This result, in particular, demonstrates the persistence of anisotropy at small scales. Like in the scalar case [ jl8| , 
dimensionless ratios of the form i?2n+i = ^n+ifi/T^^o ' 2 can be constructed using Eqs. and ([l4]) and 

evaluated at the dissipative scale t. In terms of the Peclet number Pc oc (L/i) 1 ^ we obtain 

iWocPe—, ^ +I = ^l^ti) + O(0. (15) 

The ratio for k — 2 (the so-called hyperskewness factor) diverges as Pe — > oo. This means that the conclusions drawn 



in the scalar |17| and magnetic |ll|,|l9| cases regarding the persistence of anisotropy at small scales remain valid in 
the presence of the pressure term. 

To summarize, the study of the pressure effects on the inertial-range scaling behavior of a passive vector quantity 
advected by the rapid change model reveals a double breakdown of the formal dimensional analysis. On one hand, the 
leading scaling exponents are anomalous and governed by statistical conservation laws (zero modes) of the inertial- 
range dynamics: a phenomenon already known for a local vector (magnetic) model. On the other, as a peculiar feature 
of nonlocality, the dimcnsionally regularized inertial operator develops a pole corresponding to the scaling solution 
obtained by formally matching the forcing at large scales. When the pressure term is present, nonanomalous scaling 
of the pair structure function with the dimensional exponent 2 — £ is realized only in the case when the stretching 
term in Eq. ([[]) is absent and the energy conservation enforces a constant flux solution. 
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FIG. 1. The leading exponent (upper branches) and the unphysical exponent — d + 0(£) (lower branches) vs £ for d — 3 
and A = —0.610, within the small interval A € [—0.613, —0.581] (see text). The coalescence occurs at £ c = 1.390 after which 
solutions become complex. The two branches becomes again real for £ > 1.850. 



5 



